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We use a variational method to investigate the ground-state phase diagram of a small, asymmet- 
ric Bose-Einstein condensate with respect to the dimensionless interparticle interaction strength 7 
and the applied external rotation speed Q. For a given 7, the transition lines between no-vortex 
and vortex states are shifted toward higher Q relative to those for the symmetric case. We also 
find a re-entrant behavior, where the number of vortex cores can decrease for large fl. In addition, 
stabilizing a vortex in a rotating asymmetric trap requires a minimum interaction strength. For a 
given asymmetry, the evolution of the variational parameters with increasing Q shows two different 
types of transitions (sharp or continuous) , depending on the strength of the interaction. We also in- 
vestigate transitions to states with higher vorticity; the corresponding angular momentum increases 
continuously as a function of fl. 

PACS number(s): 03.75.Fi, 67.40. Vs, 32.80.Pj 



I. INTRODUCTION 



The first experimental creation and detection of a vortex in a dilute trapped Bose-Einstein condensate relied on 
two hyperfine components in a nonrotating symmetric trap, with an applied electromagnetic field coupling the two 
components In other recent related experiments, several groups have used rotating nonaxisymmetric traps 

focussing, respectively, on the irrotational flow induced by the nonaxisymmetric shape and on the stabilization of 
vortices in such geometries |^,^,|| . Indeed, the Paris group |^,|| also measured the angular momentum of the rotating 
Bose condensate. In all these experiments, the condensate is large and generally well-described by the Thomas-Fermi 
approximation ^-jl^. 

In contrast, the present work focuses on a small rotating nonaxisymmetric condensate, which requires a small 
dimensionless interaction strength 7. Although current experiments have 7 ~ 10-1000, this value, in principle, can 
be decreased by reducing the particle number and/or tuning the scattering length with the help of a Feshbach 
resonance 13LE4j. So far, this limit has been studied for rotationally symmetric traps with a fixed angular momentum 
Lz pq-|l£ . However, since experiments do not fix the angular momentum and since we focus on nonaxisymmetric 
traps, we work at a fixed applied rotation f2. 

We first consider a noninteracting system in a rotating anisotropic harmonic trap. Although the energy eigenvalues 
for this problem have been studied previously p0[-p3| , the relevant structure of the corresponding quantum-mechanical 
eigenstates has not been considered in detail. A subset of these low-lying eigenstates of the anisotropic harmonic trap 
becomes nearly degenerate for increasing trap rotation, and conventional perturbation theory fails (as it does for a 
symmetric trap |15|| ). Therefore, we use a variational approach to study the ground state of the interacting system. 
This method permits us to determine the phase diagram with respect to the interaction strength 7 and the applied 
rotation frequency O. Monitoring the evolution of the variational ground state provides insight into the character 
of the transition to states containing one or more vortices. Moreover, we discuss the resulting angular momentum 
carried by the condensate, which increases quasicontinuously for rotations beyond that required for stabilization of a 
few vortices. 

One striking new feature of a small condensate in a rotating anisotropic trap is the possibility of increasing the 
angular momentum by elongating (which increases the moment of inertia) and simultaneously decreasing the number 
of vortices (because fewer cores can fit in the stretched form). In this way, the anisotropic condensate mimics solid- 
body rotation in a qualitatively different way from the familiar addition of more singly quantized vortices |l5[| . Our 
investigations illustrate the changing relative importance of rotational and asymmetry effects in different areas of the 
n-7 space and different values of the anisotropy parameter ujy/ujx- 

The basic formalism is presented in Sec. II along with the variational trial wave function. For a fixed anisotropy, the 
phase diagram is discussed in Sec. Ill as a function of the interaction parameter and the externally applied rotation, 
showing asymmetric vortex states and in some cases re-entrant behavior. Since angular momentum is not conserved 
in an asymmetric trap, its behavior is analyzed in Sec. IV. 
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II. BASIC FORMALISM 



In this section, we introduce the geometric configuration, construct the exact eigenstates of the asymmetric nonin- 
teracting Bose condensate in a rotating anisotropic harmonic trap, and investigate their behavior for different rotation 
speeds and asymmetries. These states provide the basis of our variational trial function for the interacting system. 
This ansatz, along with its numerical implementation, is described at the end of this section. 



A. Hamiltonian for a rotating anisotropic harmonic trap 

To study the stationary states of a low-temperature Bose-Einstein condensate in the rotating frame (when the trap 
potential becomes time- independent), we start with the time- independent Gross-Pitaevskii (GP) equation ||2^ , ^ 

+Vh- /i)* = 0. (1) 

Here iJ^*'^ = T + Vtr — is the Hamiltonian for a single particle with kinetic energy T — —h^V^/2M in an 
anisotropic harmonic trap Vtr = \M(bj'^x'^ + 'jJy'y^ + t^z^^) that rotates about the z axis with an angular velocity il, 
and Vh = sl^P is the self-consistent Hartree interaction term (the s-wave scattering length a determines the coupling 
constant g = Airah^ / M) . For definiteness, we assume tUx < tOy. 

In the absence of rotation, the noninteracting Hamiltonian iJ*^"' separates into a sum of three cartesian terms 
Hj — ^hujj{ajaj + ttjOj), where j — x, y, z. It is convenient to use the three oscillator lengths dj — ^JhjMLo] to scale 
the three cartesian coordinates, in which case the harmonic-oscillator operators have the dimcnsionlcss form 

\ ( . d \ , I 1 / d 



' 72 V dxj J V2 V dxj 

When the trap potential Vtr rotates, the term —flLz — —fl{xpy — ypx) couples and Hy, and the unperturbed 
Hamiltonian becomes = H± + Hz- Use of Eq. (||) shows that 

H±_ = ^hux (a^flx + a^al) + ^huy [a^ay + a^a^) 

+ ^^f^ [[ux + ujy) {alay - alax) + [uJx - t^y) (4^ - ayQx)] ■ (3) 

In addition to the usual "diagonal" (number-conserving) terms proportional to aja^, this operator also has an "off- 
diagonal" (number- violating) term proportional to aj,aj — ayUx- 

As in the familiar case of the two-component Bogoliubov transformation for a dilute Bose gas |2^, H± can be 
diagonalized with a generalized Bogoliubov transformation that couples all four operators. It is convenient to define 
a four-component vector a — {ax, ay, aj., ap; its elements obey the commutation relations [aj, a^] — Jjk, where J is a 
diagonal matrix with elements (1, 1,-1,-1). The Hamiltonian can now be written in matrix form as 
where 7i is a 4 x 4 hermitian matrix given by 

1 1 . (4) 
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with c = {fl/y'ujxiLJy ) I {lOx + UJy) and d = [Vl/ ^ujxUjy ) i [ujx — ujy). 

The quadratic form H± can now be diagonalized with a linear canonical transformation to a new set of four 
"quasiparticle" operators a^, defined by the matrix relation a = lAa, where the quasiparticle operators au obey the 
same boson commutation relations [aj, a[,] = Jjk ■ The transformation matrix lA follows from the eigenvalue problem 

nu^^^ ^ \kJu^^\ (5) 

where u^*^^ is the fcth eigenvector and Afc is the corresponding fcth eigenvalue obtained from the determinantal condition 
\H — AJ| = 0. The four roots are the eigenvalues [po|-p3[ 



Lol^ujl + n^T^k H - + 4c.il7^ (6) 
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where uj\ = ^(w^ + w^; these eigenvalues are identical with the classical normal-mode frequencies 

To understand our choice of notation note that the eigenvalues reduce to io\ = (wj, T for an axisymmetric 
trap (iiJx = ujy ^ In the absence of rotation (fi = 0), these two modes are degenerate. As is familiar from 

degenerate perturbation theory, applying an infinitesimal rotation breaks the degeneracy and selects out the circularly 
polarized helicity states with unnormalized wave functions ilj±{r,(j)) oc (x ± iy) exp(— ^r^) — re^"^ exp(— ir^). In 
particular, the + mode rotates in the positive sense defined by the right-hand rule and its frequency uj+ = uj± — fl 
decreases with increasing angular velocity as is obvious when viewed in the rotating frame 

In the general case of an anisotropic trap with Ux < Wy, the two modes lu± are nondegenerate even for $1 = 0, 
when the + mode with uj-^ — lo^ is linearly polarized along x and the — mode with w_ — LOy is linearly polarized 
along y. For nonzero rotation speed fJ, the dispersion relation (|^) exhibits a crossover near ^ [lo^ — lo'^) / Auj ± from 
an anisotropy-dominated regime that favors linearly polarized states at small to a rotation-dominated regime that 
favors elliptically polarized states with definite helicity at large fl. After the cross-over (at relatively small rotation 
speeds, for example J7 « 0.05 oj^ for an asymmetry of 10%), the smaller positive eigenvalue ujj^ decreases nearly linearly 
for moderate fl; in the large-0 regime (fl/ujx ^ 1), however, vanishes like {ujy — u!x){(^x — ^) with an infinite 
slope a.s ft —^ ujx- In contrast, the larger positive eigenvalue lj_ increases with increasing Q,. Figure ^ illustrates the 
positive eigenvalues uj± as functions oiQ for two asymmetries. For the small asymmetry LOy = 1.014 Wj^, the deviations 
of u!± from the axisymmetric case are hardly visible. For ujy = 1.1 cj^:, the deviation of cOx from a straight line is seen 
at both ends of the range of rotation speed < fl < uix- 

The four eigenvalues in Eq. (^) of H± in Eq. (||) can be taken to constitute a diagonal matrix A with elements 
(tj+jCij-, — —L0-). For the two positive eigenvalues uj±, the eigenvectors have the form 



/ cosh x± cos 9± 
±i coslixi sin 9± 
— sinh x± sin ri± 

\ ±i sinh x± cos r?± , 



(7) 



These two eigenvectors satisfy the normalization condition u^-'-'t Ju*^*^) — Sjk, where + and — correspond to j and 
k — 1 and 2, respectively. It is not difficult to obtain explicit expressions for the (real) hyperbolic parameters x± 
and for the (real) trigonometric parameters 6± and 77±; they depend on the trap frequencies ujx and Uy and on the 
external rotation speed fl. Symmetry considerations readily relate the remaining two eigenvectors (those for the two 
negative eigenvalues — w±) to u^^\ The resulting four eigenvectors u^^' with j ~ 1, • • • ,4 form a complete basis set 
and obey the normalization condition 

yO-)t Jy(fe) ^ J^^ . (8) 

(k) 

The transformation to the quasiparticle operators is determined by the matrix Ujk = Uj of the four eigenvectors 
written in successive columns (it is the analog of the "modal matrix" that plays a central role in the theory of 
small oscillations of mechanical systems about stationary configurations [^). In this way, the matrix U satisfies the 
eigenvalue equation [compare Eq. (^)] 

nU = JUA, (9) 

and the normalization (^ for the eigenvectors implies that JU = J. Equivalently, multiplication of the eigenvalue 
equation (|p) by gives U^TiU = JA, showing that the transformation matrix U indeed diagonalizes the hamiltonian 
matrix Ti. with the appropriate boson metric J. Correspondingly, the quasiparticle operators follow from the matrix 
equation 

a = JU'fJa, (10) 
and it is easy to see that the transformed unperturbed Hamiltonian has the expected diagonal form 

H± = ^ha^Tia = ^hijj^{a\_a+ + a+a^) + ^huJ-{aLct- + a^a^___) (11) 
that represents a set of uncoupled harmonic oscillators. Evidently, the spectrum of allowed states has the eigenvalues 

fn+n- = + |)^+ + Hn- + ■^)^-, (12) 

where n± is a nonnegative integer. As is clear from Fig. |^, the set of lowest eigenvalues e„^ = h{n^ + ■^)w+ vanishes 
as n approaches the confinement limit (fl lox). This degeneracy precludes a simple perturbation approach for the 
interacting system. 
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Equation ( |To| ) yields explicit expressions for the quasiparticle operators 

a± = coshx± (cos 9±ax T 'ism9±ay) + shThx± isi'nij±al. ± icosri±a^y) , (13) 

along with their adjoints. Here, the hyperbolic parameters x± the trigonometric parameters 6± and ri± guarantee 
the correct commutation relations for the quasiparticle operators. In particular, the parameters x± determine the 
number- violating "Bogoliubov" coupling between the a and operators. As expected from the form of Eq. (|^), 
these parameters vanish for = 0, and x- remains small and positive for all allowed > 0, whereas x+ is negative 
and decreases rapidly as lOx- For an axisymmetric trap, x± vanish identically for all fl, and the quasiparticle 
operators a± then reduce to the familiar helicity operators a± = {ax T 'icLy)/\/2. 

We can now analyze the noninteracting eigenstates. The rotating ground state ipoo{x,y) is determined from the 
pair of conditions a±ipQo — 0; it has the form 

foo{x, 2/) = ( j exp [-^{Ix'^ + 2imxy + ny^)] , (14) 

where the real quantities I, m, and n depend on the trap frequencies and the rotation speed |2l| ]. For an axisymmetric 
trap with uox = = the ground state has the expected isotropic structure with I = n — 1 and m — 0. This state 
remains isotropic for all Q < luj_. 

For an anisotropic trap with ujx < '^yi however, the noninteracting ground state has a nontrivial phase oc xy that 
represents the irrotational flow induced by the rotating trap ||2^ , |26| , |l0[p^ ,^||30|] . For small rotation speeds and small 
asymmetry, the ground-state density is indistinguishable from that of the symmetric trap, so that it has an essentially 
circular shape fitting wholly into the elliptical trap geometry. For fast rotations, when the angular momentum per 
particle is considerable, the condensate experiences a torque that stretches the ground state density along the axis of 
lesser confinement (the x axis in our choice). For 10% trap asymmetry, the central peak also decreases appreciably. 
As VL approaches the smaller transverse oscillator frequency (here, ujx)^ the parameter / vanishes. Thus the oscillator 
confinement for the noninteracting anisotropic trap disappears at the critical value 17 = Wa;, as is well known from 
classical mechanics. 

The normalized excited states are given by the familiar harmonic-oscillator construction 

{alr+ {air- , , 

^n+n. (x, y) = — == -j==- <Poo(a;, y)- (15) 



In the axisymmetric, static limit {ujx = Wj,, f2 = 0), tpio and (poi reduce to the degenerate pair i}}± discussed below 
Eq. (^). For ujx < i^y, the first excited noninteracting state tpiQ has an excitation energy and represents a 
nonaxisymmetric vortex with unit positive circulation and a node at the trap center. The other noninteracting singly 
quantized vortex state (p^i also has a node at the trap center, with higher excitation energy hio^ and unit negative 
circulation. 

A recent study of a small rotating axisymmetric Bose-Einstein condensate works at fixed angular momentum 
Lz = hi (which is appropriate only for an axisymmetric trap). The resulting equilibrium configuration is then 
determined by minimizing the total energy iV£'iab(0 subject to the constraint of fixed I. The corresponding angular 
velocity $1 then follows from the relation 51 — dE^h/hdl. 

The constraint of fixed = Kl is analogous to the constraint of fixed total N in the canonical ensemble. As in the 
transition to the grand canonical ensemble, however, it is often advantageous to eliminate the constraint by making 
a Legendre transformation from fixed / to fixed 57, which here merely means transforming to a rotating frame. Since 
the resulting Hamiltonian H = H^^^^ + Vh then contains the term —V.Lz, the expectation value (L^) for the angular 



momentum as a function of follows directly from the Hellmann-Feynman theorem |31| (L^) = —dE{n)/dVt, where 
E = (H) is the energy in the rotating frame. For example, the expectation value of the angular momentum for the 
noninteracting eigenstate (p„^„_ is simply 

Ln+n^/h = -{n+ + \)duj+/dn - (n_ + \)dw-/dn, (16) 

as follows directly from Eq. ([T^). Figure |l] shows the dependence of u}± on the external rotation U,, and the resulting 
Ln+n^ also depends on Q for any nonzero trap anisotropy. With the definition L± = —hduj±/dn, we have Lqo = 
i(Z/+ + L_) and Lio = Lqq + L+. These angular momenta per particle for the two lowest noninteracting states (poo 
and ipiQ are included in Figure For both anisotropics shown {cuy/ojx = 1.014 and 1.1), Loo remains small until 
n approaches lOx', in contrast, Lio rises rapidly and linearly for small fl, remains close to one quantum of angular 
momentum for most of the allowed range, and then grows rapidly as Ux- 
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The asymmetric noninteracting excited states possess a rich structure and the complete density distribution is 
needed to characterize them. As an example, Fig. y presents density contours of |<P3oP across the whole xy plane, for 
slow, medium and fast rotation and for two different values of ujy/ojx > 1. The asymmetry aligns the vortex cores 
along the axis of weak confinement. At slow rotation {^l/uj^ — 0.01), most of the condensate density accumulates 
at the ends on the weak axis. For the case of ojy/uJx — 1-1, the depletion along the tight axis essentially splits the 
condensate at the position of each vortex "core." This result is plausible because this small fl is in the asymmetry- 
dominated regime where fluix <C ujy — w^. For medium rotation {Q,/ujx = 0.5), the effect of the asymmetry smears out, 
and the densities resemble rotation-distorted symmetric trap eigenstates, with the vortex cores located close to the 
trap center and surrounded by a region of nearly constant density toward the edge. For fast rotation [VI/lOx — 0.95), 
the condensate stretches along the weak axis, particularly pronounced for ujy/oJx = 1-1. There, the non-ellipsoidal 
shape of the inner most contour indicates that two of the three overlapping vortex cores are now off-center on the 
X-axis. The corresponding loss in angular momentum is again compensated by having a larger density at the ends of 
the weak axis. Note that the depletion along the tight axis is here much less than for slow rotation so that the core 
area is still fully surrounded by the condensate. 



B. Variational ansatz 



The noninteracting eigenstates now allow us to treat the interacting problem. With the previous dimensionless 
spatial variables, the GP equation (|l|) becomes 

{H^^^ - H + ^'K-ih,/U^^l^*il))ij = Q, (17) 

where — H± + Hz, with H± in diagonal form taken from Eq. (pl|), 7 = Na/dz is the small interaction parameter, 
and "0 is the condensate wave function normalized to 1. Note that we assume a nonaxisymmetric trap rotating 
with a fixed angular speed f2. This approach describes a "helium-bucket-like" experiment pl o] and complements 
the theoretical work on vortices in small axisymmetric condensates (for example, Refs. Jl5| , p^ ), where the angular 
momentum (a good quantum number) is fixed. 

As is obvious from the vanishing of the noninteracting eigenfrequency w_|_(il) for ^ ^ uix, the lowest eigenstates Lpno 
become nearly degenerate for asymmetric traps, and straightforward perturbation theory is not applicable. Because 
we are interested in small condensates with 7^1, however, the noninteracting eigenfunctions discussed in Sec. II. A 
provide a suitable basis set for an expansion of the planar part of the interacting condensate wave function. Our 



strategy is thus to use a linear combination of the lowest eigenstates as a trial function |15|, 



il)n{x,y) =^Csipso{x,y) , (18) 

s=0 

where n indicates the cutoff of the nth-order trial function at the excited state </3„o- This cutoff makes the calculation 
tractable; it can be chosen so that the trial function captures the essential physics for rotation speeds at which higher 
excited states do not contribute. This trial function restricts the region in the phase diagram that we can investigate. 

We assume that the bosons are in their ground state along the axis of rotation (po(^) — 7r~^/* exp(— 2^/2). The 
variational ground state follows by minimizing the free-energy functional 

= ^ IcsP^i"'' + V2Ti-f^ujxUJy c*c*CkCiIijki , (19) 

s ijkl 

with respect to the set of variational parameters . Herews*^' ~ (s + ^)aj+ + ^(ci;_-|-Wz) and Zy/i;; — J d'^rip*Qip*Qipf^fjipiQ. 
The variation is constrained by the normalization condition 

/„ n n 

d^rip*il}= I d^r ^ clc^,iplf^{x,y)ip,,Q{x,y) ^^\cs\^ ^l. (20) 

ss'=0 s=0 

In general, the parameters Cs are complex, leading to 2 x (n-|- 1) real variational parameters. This number is reduced 
by one through the norm condition ( |20| ) . We can eliminate another degree of freedom by fixing the arbitrary phase 
of the wave function. Therefore we have to minimize the energy ([l9| ) for an nth-order trial function with respect to 
2n independent variables. 
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C. Numerical implementation 



We use a well-known simplex algorithm [ p2[ to minimize the energy ( |l9| ) in terms of the Cs as variational variables. 
The algorithm basically determines the function to be minimized with respect to 2n variational parameters at the 
corners of a 2n + 1-dimensional simplex in parameter space. We reduce the number of independent variables by 
choosing the phase of one of the Cs to be equal to zero and rewriting its modulus explicitly as a function of the other 

variational parameters |cfe| = ^/l — J2s^k I'^sP- "^^^ energy functional is restricted to the unit sphere. We filled the 
space outside by taking the values of the energy on the surface of the unit sphere and magnifying them with increasing 
distance from the sphere. This procedure allows us to minimize in an unrestricted space and yet ensures that the 
variational equilibrium state satisfies the norm condition. In principle, it makes no difference which Ck we exclude 
from the minimization procedure. However, since we have to find the minimum with respect to many parameters and 
since the energy functional has many local minima, we swept all fc = 0, . . . , n for a given point in the phase diagram. 
This way, we repeated the minimization on n + 1 different representations of the energy and thus drastically enhanced 
the likelihood of finding the true global minimum. 

We monitored the reliability of our results by comparing the contributions of higher excited states. We conclude 
that we have captured all important ingredients when the distribution of weights among the Cs was not affected by 
including higher excited states into the trial function. For reasons of symmetry, this test required increasing n by 
at least 2. The minimization was very robust and reliable for most rotation speeds and interaction strengths that 
we examined. The behavior of the resulting Cs with increasing rotation f2 was generally very smooth apart from the 
transition lines between equilibrium states with different numbers of vortices. Another check was the convergence of 
the equilibrium energies for different n. If a higher basis state was important, the energy clearly decreased when it 
was included. If the additional dimensions in the minimization space were irrelevant, the resulting minimal energy 
was the same within the numerical errors. 



III. THE PHASE DIAGRAM 



In this section, we determine the lowest transition lines between variational equilibrium states in Q-j space for 
symmetric traps and for weakly and moderately asymmetric traps. In discussing the variational equilibrium conden- 
sate, we clarify our criterion for the transition with increasing rotation speed il when the interaction and the trap 
geometry are fixed. Moreover, we investigate the change in the character of the transition in different regions of the 
phase diagram. 



A. Symmetric geometry 

First, we review the phase diagram for an axisymmetric trap that serves as a comparison for the new features induced 
by the asymmetry. The transition lines between equilibrium states are determined by minimizing the energy ( p^ ) for 
many sets of parameters {ujy,fl,"f) and comparing the resulting states. We call a state a q vortex if ipqo dominates 
the complete trial function ip, namely |cgp > |csp, for all s ^ q. A transition between two states is identified by a 
change in the dominant weight from |cgp to |c<j'p. This criterion works well for "clear" states, where one particular 
|cqp dominates strongly over the other contributions. 



From previous theoretical studies [|5 33 1 and from the experiments of the Paris and MIT groups [|6|,[7|, we expect 
that the axisymmetric condensate accommodates the angular momentum associated with the rotation by a sequence 
of transitions to states with more vortices. For the axisymmetric trap potential, our results agree completely with 
the earlier theoretical phase diagram . As seen in the dot-dashed lines in Fig. |^, the first two transition lines are 
strictly linear in 7 and pass through Q — ujx io. the noninteracting limit (7^0). As expected, the line f2c(7) for the 
first transition to a state with a single vortex agrees exactly with the critical frequency i^c/i^x = 1 — "f/VSn obtained 
with first-order perturbation theory [ p"5|j3^ . Most of our data are obtained with an 8th-order trial function (the area 
very close to uj^ has been tested by using n = 10). For the phases with up to two vortices, only the lowest five 
iPsQ contribute significantly. 

In order to characterize the various phases more precisely, it is instructive to monitor the behavior of the variational 
parameters along a vertical cut in the phase diagram, namely with increasing Q for a fixed interaction strength 7. Such 
a cut is included in Fig. ^for an axisymmetric trap with 7 = 1 (we measure £7 in units of uJx)- We plot the |csp starting 
in the no- vortex phase. Below the range shown in Fig. H, we always find |coP = 1, confirming that the noninteracting 
ground state is also the variational equilibrium state for these slow rotations. At flc = 0.80053 = 1 — I/VStt, there 
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is a sharp transition to a new state that consists purely of the noninteracting ipiQ state. This behavior identifies the 
critical frequency for thermodynamic stability of a singly quantized central vortex. The next three transitions to 
other combinations are also clearly seen, indicating a sequence of transitions to states with an increasing number of 
singly quantized vortices, each in a well-separated range of Q. Beyond the one-vortex phase, the variational states 
are mixtures of various noninteracting states. In particular, the two-vortex phase involves mixing with other states 
of two- fold symmetry (v3oo,¥'4o) and the three- vortex phase similarly contains other states with three-fold symmetry 
{(poo, (peo)- For rotation speeds higher than « 0.98 Wa;, we are too close to the degeneracy limit D, — lox to exclude the 
possibility that higher-order trial functions might change the distribution of weights. In the present symmetric case 
(dot-dashed lines in Fig. 3), the order of appearance of the different phases is the same for any vertical cut at fixed 7 
in the phase diagram. 

We have included the angular momentum L^i^l) = —dE/dfl in Fig. ^ for 7 = 1. Each transition that adds one 
more vortex induces a discontinuous upward jump in the angular momentum. Evidently, not all values of angular 
momentum are allowed. Specifically, values in the range < < 1 are absent, which is consistent with the 
character of the equilibrium functions before and after the first transition (they consists purely of foo or 9310, carrying 
exactly zero or one quantum of angular momentum respectively). For the allowed ranges of Lz, direct comparison 
with Ref. jig]] shows that we found the same mixtures of (pso, although they fix the angular momentum. Since the 



transition lines are linear in 7, we can also make contact with the results from Ref. |15 . In particular, the first three 



transition frequencies indicated in their Fig. 2 for fixed angular momentum are the same as those found here for fixed 
rotation speed f2. Furthermore, the phases have the same symmetry, the lowest two being pure ipQQ and ipio and 
thus rotationally symmetric, whereas the equilibrium state of the q ^ 2 and q = 3 phases have off-center vortex cores 
arranged to give a two-fold or three-fold symmetric structure. 



B. Asymmetric geometry 

We investigate two specific trap asymmetries in detail. A weakly distorted trap with ujy/cux = 1.014 already shows 
some new features and reflects a delicate balance between the symmetric (rotation-dominated) and asymmetric (trap- 
dominated) influences. A second trap geometry of ujy/ujx — 1.1 displays more pronounced effects of the asymmetry 
and thus provides a clearer picture of the basic physics. Figure |^ shows the first two transitions in the j-n plane for 
both asymmetries, along with the corresponding curves for the symmetric trap with ujy/uj^ = 1. 



1. Transition lines 



First, note that the transitions for the asymmetric traps occur at higher rotation speeds U, than in the symmetric 
case; in addition, the shift increases with increasing asymmetry. In contrast to the linear behavior (fi^ oc 7) of the 
symmetric system, the transition lines curve significantly for small interaction strengths and high rotation speeds. 
Most remarkably, there is a critical threshold coupling constant 7c below which a ground state with a singly quantized 
vortex is never favorable. This behavior is understandable because the vortex core size decreases with increasing 
interaction parameter and only fits into the trapped condensate for not too weak interactions. Since the semi-minor 
axis of the ellipsoidal trap fixes the size available for a vortex core, a smaller asymmetry allows the introduction of a 
vortex at lower 7. 

Second, the term —VLLz in the Hamiltonian tends to favor states with large angular momentum. In addition to 
increasing the number of vortices, a greater trap asymmetry makes the condensate more susceptible to rotation- 
induced elongation, placing a greater part of the condensate farther away from the rotation axis, increasing the 
moment of inertia and hence the angular momentum. In this way, the condensate can accommodate a higher angular 
momentum without introducing vortices. 

For fixed interaction strength 7 not too far above 7c, we even find a re-entrant region where the one- vortex phase 
is followed for higher rotation speeds by a no-vortex phase. The tip of the one-vortex phase surrounded by the 
no-vortex phase for the small asymmetry [ujy/uix = 1.014) is illustrated by the occupancies |csp along a cut through 
that tip (at 7 — 0.08935, Fig. ||). The occupancies in the re-entrant no- vortex phase continue as if it had never been 
interrupted by the one-vortex phase. For illustration. Fig. 5 also includes typical density contours for J7 around the 
lobe tip. The elongation along the horizontal x axis for the vortex-free states is pronounced. The thinner waist of 
the condensate on the vertical axis results from the admixture of ip2o- The one- vortex state again has only Lpio as a 
constituent and illustrates how the circular ring of maximal density in the symmetric case deforms to two pronounced 
density peaks on the x axis for this range of fi. This re-entrant behavior reflects the singular character of the limit 
f2 — > oj^ for asymmetric geometries, when the confinement parameter / in Eq. (q3) tends to zero. Below the threshold 
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interaction strength for vortex stabilization, the density contours always represent elongated no-vortex states; their 
width is smaller for smaller interactions strengths and increases slowly with increased rotation. This elongation seems 
to hinder vortex formation. 

In order to study the detailed structure of the phases in the asymmetric cases, we considered again the intermediate 
interaction strength j — 1 and determined the occupancies |csp for small {ujy/ujx = 1.014) and moderate {uiy/ujx — 1.1) 
asymmetry, as shown in Fig. ^jand Fig. ^ respectively. As the applied rotation increases, the smaller asymmetry shows 
phases with increasing number of vortices, just as in the case of a symmetric trap (compare Fig. ^ for a symmetric 
trap). In contrast to the re-entrant behavior for 7 = 0.09 (Fig. |^), we conclude that these relatively strong interactions 
7 = 1 eliminate the effect of the (small) asymmetry, in part because the transitions occur at slower angular velocity. 
For the moderate asymmetry [ujy = l.lujx), however, the phase diagram in Fig. |^ still exhibits re-entrant behavior; 
thus a 10%-asymmetry dominates the behavior for this interaction strength (7 = 1) and precludes more than two 
vortex cores. Note that the one- vortex phase continues for fl/uj^ > 0.99 as if it had not been interrupted by the 
two- vortex phase. 

Typical density contours for the two asymmetries are included in Fig. ^j. Several features differ significantly from 
the symmetric geometry. The no-vortex phase has a considerable ip2o admixture, producing a constriction along 
the y axis. In the one-vortex phase, the condensate has a central vortex, but the elliptical trap and rotation-induced 
elongation deform the condensate noticeably. For LUy — 1.1 uj^^ the contributions from ip^Q (and other odd states) grow 
as n approaches the transition to two vortices, deforming the surface region because of four vortex cores that move 
in from infinity. The density contours for the interacting two-vortex state show the admixture of the noninteracting 
ground state because the vortex cores are pushed further from the center of the trap than in the symmetric case. The 
small admixtures of other basis states with an even number of vortex cores favors the accumulation of density closer 
to the center (outside the core regions). 

For the smaller asymmetry, the subsequent figures display three and four separate cores, respectively. In the three- 
vortex phase, however, the cores are not symmetrically distributed around the center (the appreciable occupation 
of (^20 enhances the two vortices along the x axis, placing the third core on the horizontal axis and further away 
from the trap center). From Fig. ^, note that the mixture of noninteracting states differs from the symmetric case 
(Fig. ^), where the states contained only noninteracting states with the same rotational symmetry. The last two 
density contours in the right column illustrate that the re-entrant phases can indeed be classified as states with one 
or zero vortices. For such fast rotations, the condensate is very elongated and very fiat (as can be seen from the fewer 
density contours) . It is energetically favorable to reduce the number of vortices (which need a wider condensate) and 
to compensate the loss of vorticity in the vortex cores by expanding the condensate along the horizontal axis. In the 
re-entrant ipoo-pliase, we have several small contributions from basis states with an even number of cores. This causes 
small ripples on a thin extended Gaussian density. Basis states with an odd number of cores would put a density 
minimum at the center of the trap, but the dominant ipoo suppresses this tendency. 

In recent experiments the Paris group measured the nucleation of vortices in their large, cigar-shaped condensate 
and the corresponding angular momentum resulting from the vortices alone. When ramping up the rotation, they 
eventually find vortex-free states again, which is a re-entrant phenomenon similar to what we find in our analysis for 
much smaller condensates. Moreover, they also measure the nucleation of vortices when sweeping the asymmetry e 
and keeping the rotation fixed, leading again to a window of vortex stabilization. For the small condensate, we also 
determined the transition to a one-vortex state in the Cl-LUy plane, as illustrated for several 7 and asymmetries up to 
20% ( Fig. ^). For asymmetries larger than the rightmost end of the graphs, there is no vortex state for rotations 
up to ^/tOx = 1. In agreement with the Q-^ phase diagram (Fig. ^), we here see the re-entrance in that Q-ujy phase 
diagram, leading to a qualitatively similar restricted rotation window for stabilizing a vortex. 

2. Types of transition 

Having discussed the typical phases, we can now consider the details of the transitions themselves. The occupations 
of the separate constituent noninteracting states (compare Fig. ^ and Fig. ^ have an important new feature. There 
is a smooth transition between the no-vortex and the one-vortex state for 7 = 1 for both asymmetric geometries, 
reminiscent of a second-order transition. This phenomenon can already be found for a very small asymmetry, as 
illustrated for ujy = 1.001 w^; in Fig. ^ The critical frequency Qc = (0.80125 ± 0.00005) w^; is a few per cent larger 
than flc = 0.80053 for the symmetric trap. Nevertheless, the density contours close to the transition reveal that a 
vortex core gradually enters the condensate along the y axis (as seen in Fig. [l0| ). This behavior is qualitatively distinct 
from that for the symmetric case (for the first transition at 7 = 1, we examined points as close as rotation speeds 
A{n/ujx) — ±10^^). Thus we infer that the character of the transition in asymmetric traps differs fundamentally 
from the symmetric case, for we now have a cross-over region where both cq and ci are nonzero. In this situation. 
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the meaning of the critical transition frequency fie becomes somewhat blurred. For both asymmetric geometries 
(ujy/ujx = 1.014 and 1.1), this cross-over region shrinks for smaller interactions. In fact, we find a sharp transition 
below 7 < 0.1 for the smaller asymmetry and below 7 < 0.8 for the moderate asymmetry. For all other subsequent 
equilibrium phases, we found spontaneous jumps in the occupancies, similar to first-order transitions. 

The change in character can be understood as follows: for parameters that lead to a sharp transition, the energy 
functional has two competing, well-separated main minima, one lying in the co-dominated sector and the other in 
the ci-dominated sector. Indeed, depending on whether we choose cq or ci to implement the norm condition in our 
minimization procedure, we find one or the other minimum in the neighborhood of the transition. The comparison 
of the energies then gives the true global minimum. The depth of these two minima gradually changes with fi; at 
the transition, the minimum representing the one-vortex state becomes deeper. This picture allows for hysteresis in 
stabilizing a single vortex with increasing rotation. In fact, hysteresis is the favored explanation for the deviation 
of the measured flc from the Thomas-Fermi predictions (see p|,p|,p5 36 1 and references therein). Moreover, the fact 
that the phase after re-entrance appear to be the continuation of the phase before the previous transition means that 
the corresponding minimum still exists and again lies below the energy of the intervening higher-vortex state. For 
a continuous transition, in contrast, the minimum energy functional must lie in a valley connecting the two sectors. 
Here, the different implementations of the norm condition lead to the same minimal state. With increasing Q, this 
global minimum gradually moves along this valley from the co-dominated sector through a cross-over region to the 
ci-dominated sector. 

From Fig. |^ for the smaller asymmetry with ujy — 1.014^3;, we can also observe that the three-vortex phase 
undergoes much stronger changes across its range than other phases. Although there is a sharp transition to a phase 
with (/930, an appreciable amount of the two- vortex basis-state remains and 1^920 dies out only gradually. This behavior 
yields a very complicated picture for smaller interaction strengths (roughly at 7 « 0.6), where the </33o-dominated 
phase disappears completely leading to a direct transition from two to four vortices. For even smaller 7, the two- vortex 
phase is followed by a one- vortex phase before the four-vortex phase develops. In part, this complicated picture arises 
from the suppression of the 1^930 contribution to the one- vortex phase (which in the moderately asymmetric trap with 
ujy = l.lujx leads to an contribution of 9550 rather than 1^930 to the one- vortex phase). 



IV. THE ANGULAR MOMENTUM 



The asymmetry in the trap geometry breaks the cylindrical symmetry, so that the angular momentum Lz around 
the axis of rotation is no longer a good quantum number. We therefore investigate the effect of the asymmetry on Lz 
in some detail. 

The angular momentum (in units of h) as a function of the trap rotation fl is shown in Fig. ^for all three geometries 
{ujy/ujx = 1.0, 1.014, and 1.1) for the fixed interaction 7=1. The intimate relation to the occupancies within the 
various phases for each geometry can be seen in the corresponding plots (Figs. |^, ^, and 0). The sharp kinks in Lz 
for the symmetric trap arise from the sudden changes in the occupancies and the different angular momentum carried 
by the noninteracting states. The finite slope of the plateaus beyond the second transition refiects the mixing of the 
various noninteracting states in the variational ground state and the off-center positions of the vortex cores |1^ . 

For the asymmetric rotating trap, even the noninteracting eigenstates ipso with a central s-fold vortex carry an 
angular momentum different from s (in fact, the angular momentum diverges for — > lOx). The divergence comes 
from the contribution of the circulating quanta with positive helicity, i+(r2), occurring in the angular momentum 
of every basis state ipsO- This is easily seen by expanding = ^duj^/dil for small lUx ^ ^ and small asymmetries 
= 1 -I- 2e^, leading to 

L+«^=L= torn^iOx. (21) 

y/LJx — il 

In asymmetric traps, this divergence dominates the angular momentum in all noninteracting eigenstates ipso for high 
rotation speeds and hence in any linear combinations of them. Thus, even discrete changes in the structure of the 
equilibrium state (as at a transition) will have a far less dramatic effect on the angular momentum at large rotation 
speeds near tOx- Moreover, asymmetric condensates are significantly elongated along the axis of weaker confinement 
as f2 — !■ Wj.. The resulting redistribution of density (and mass) with respect to the axis of rotation produces an 
additional, nonquantized contribution to the angular momentum. 

Compared to the symmetric condensate, the angular momentum for the trap with smaller asymmetry at 7 = 1 
has its jumps smoothed over a small but finite range of Q because of the continuous transition where one off-center 
vortex gradually moves towards the center (see Fig. in addition, the transition is shifted toward higher rotation 



speeds (see Fig. 11). The subsequent sequence of phases corresponds to an increasing number of vortex cores that 
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cause nearly vertical jumps just as in the symmetric case. The divergence of the angular momentum dominates only 
at ^/lOx ~ 1- 

For the moderate asymmetry, the stronger admixture of ip2o in the no- vortex phase yields a faster rise of the angular 
momentum below the first transition. The first kink is even more tilted because of the extended range for the transition 
to a one- vortex state. Inside this one- vortex phase, there is a broad region with a pure ipiQ state (up to fi/ciJa; ~ 0.92). 
The angular momentum, however, exceeds one because of both the initial growth of with 17 and the deformation 
of the condensate density. The rapid rise of Lz below the second transition reflects the growing admixture of ip5Q. 
That second transition, seen as a discontinuous derivative in Lz, occurs at an angular momentum Lz ~ 2.4, already 
exceeding the value for a symmetric trap with two vortices. Beyond that point, the angular momentum increases 
smoothly with fi, even beyond the re-entrant transitions at « 0.9905^3; and ~ 0.9945 respectively (see inset 
in Fig. pi]). In this limit, the divergent completely dominates the behavior. 



V. CONCLUSION 



The behavior of a small Bose condensate in a rotating anisotropic trap differs significantly from that of a symmetric 
condensate, even for very small asymmetries. In large part, this difference arises from the stretching of the condensate 
along the axis of weak confinement, especially for lo^- In addition, it reflects the irrotational flow induced by the 
rotating confining potential that pushes on the gas when viewed from the laboratory frame |p9|,p 

For any specific asymmetry and sufficiently small interaction strength, a one-vortex state can never be stabilized, 
as illustrated in Fig. y. For somewhat larger values of the coupling constant 7, the first transition, from a no-vortex 
state to a one-vortex state, is sharp, but the subsequent transitions are re-entrant, back to a no-vortex state with 
significant elongation that carries the relevant angular momentum (see Figs. || and ^). For these weak interactions, 
the transitions all occur close to « w^:, so that the stretching dominates. 

For larger interaction strength (7 ^ 1), the transitions occur at lower values of fi, and they become quasicontinuous 
in that the occupation |cop of the noninteracting ground state vanishes smoothly (see Figs. ||, 0, and [l0| ). As in the 
case of a classical fluid in a rotating elliptic cylinder p9| , the angular momentum is finite even below the transition 
to a state with one vortex because of the irrotational flow induced by the rotating asymmetric trap (see Fig. |l^) . 
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FIG. 1. The positive eigenvalues lj± and the angular momenta Loo and Lio of the two lowest noninteracting eigenstates 
tpoo and tpio as a function of the rotation speed Q, for asymmetries Uy — 1.014 (dashed lines) and cuy — 1.1 (solid lines). All 
quantities are in dimensionless units (scaled with respect to the lower oscillator frequency ujx)- 



FIG. 2. Density contours of the noninteracting |(^3o|^ across the xy plane for the traps with ujy/ujx = 1.014 (left) and 
ujy/cux = 1.1 (right) for Q/u^ = 0.01, 0.5 and 0.95 (top to bottom). Distances are scaled in units of dx in both directions and 
the width shown is 6dx across. 



FIG. 3. The two lowest transition lines for uiy = 1.0 uJx (dot-dashed lines), ujy = 1.014 cj^ (dotted lines) and uiy = 1.1 0;^; 
(solid lines), determined with trial functions up to 10th order. The lowest transition lines (a,b,c) represent the critical rotation 
Qc for the stabilization of a single vortex; at the second transition lines (d,e,f) the condensate starts to be dominated by the 
(^20 ■ There are unresolved higher vortex-phases beyond the second transition for sufficiently large interaction strength. 



FIG. 4. ujy — 1.0 ujx, 7=1: Occupancies (left scale) and (right scale) for axisymmetric condensate as function of fl for 
fixed interaction 7 = 1. The numbers above the lines denote the corresponding |cs]^ and the starred solid line represents the 
angular momentum L^. 
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FIG. 5. Occupancies |csp for the basis-states for iOy — 1.014 oj^,, 7 = 0.08935, calculated with a 8th-order trial function and 
typical density contours (at n = 0.99515, 0.9952, 0.9953). Contributions for s > 2 are negligible. The tip of the pure ipio-lohe 
cuts sharply into the no-vortex phase, which has an appreciable 1/P20 admixture. Note the narrow range in SI. The density 
contours are shown over 6 oscillator lengths dx in each direction. 

FIG. 6. ujy — l.OlAuix, 7=1: Occupancies (left scale) and (right scale, note upward shift). The numbers above the lines 
are short for the correspondent |csp and small contributions from \c(i\^ . . . |ciop are not separately labeled. The starred solid 
line represents the angular momentum L^. 

FIG. 7. Uy = 1.1, 7 = 1: Occupancies (left scale) and Lz (right scale). The numbers above the lines are short for 
the correspondent \csf and small contributions from jcgp, |c6p...jciop are not separately labeled. The starred solid line 
represents the angular momentum L^. 



FIG. 8. Density contours for variational equilibrium states for small (1.014 left) and moderate (1.1 right) asymmetry and 
7 = 1. The pictures represent typical states within the different phases (cf. Fig. ^ and Irrespectively). For the small asymmetry, 
the number of vortices increases from zero to four (n = 0.805, 0.935, 0.95, 0.97, 0.982) whereas in the moderate asymmetry, 
we find re-entrant behavior into one and zero- vortex states (Jl = 0.87, 0.9, 0.97, 0.993, 0.997). The maximum width shown is 
6 oscillator lengths d^. 

FIG. 9. Lines of critical rotation for the stabilization of one vortex versus the asymmetry, for 7 = 0.8, 1, 1.5 (left to right), 
from 8th-order trial functions. The one-vortex states are to the left of each graph, the vortex-free ones to the right. All 
quantities are given in dimensionless units. 



FIG. 10. Occupancies and density contours for ujy = 1.001, 7 = 1, from an 8th-order trial function. The occupancies for ipso 
and higher are not labeled for clarity. The density contours are taken at rotations Q. = 0.8007, 0.8011, 0.8012, 0.8013 and are 
about 4 oscillator lengths d^ across. 



FIG. 11. The angular momentum Lz in units of ft as a function of the angular velocity Q for the three geometries 
ujy/ujx = 1.0, 1.014, 1.1 (dot-dashed, dotted and solid lines respectively). The inset shows the magnified re-entrance range 
for uiy/uix = 1.1. 
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